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Abstract 

We study the transport and locaHzation properties of scalar vibra- 
tions on a lattice with random bond strength by means of the transfer 
matrix method. This model has been recently suggested as a means 
to investigate the vibrations and heat conduction properties of struc- 
tural glasses. In three dimensions we find a very rich phase diagram. 
The delocalization transition is split, so that between the localized 
and diffusive phases which have been identified in the Anderson prob- 
lem, we observe a phase with anomalous, sub-exponential localization. 
For low frequencies, we find a strongly conducting phase with ballistic 
and super-diff'usive transport, refiecting a diverging diffusivity. The 
last phase generates an anomalous heat conductivity which grows with 
the system size. These phases are the counterparts of those identified 
in an earlier study of the normal modes. 



1 Introduction 



It is well known that disorder in the parameters of physical models can 
induce scattering and localization, and the most outstanding example is the 
electronic Anderson model. This model has been the focus of a large set of 
theoretical, numerical and experimental studies over the last three decades, 
and there is an enormous literature on various aspects of this problem. At the 
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same time, the related problem of vibrations in disordered materials received 
almost no attention, although structural disorder is very common in glasses 
and other disordered materials. 

It is well known that in low temperatures there are several features com- 
mon to most glasses. An outstanding phenomenon is the existence of an 
anomalous heat conductivity at low temperature, a plateau in the heat con- 
ductivity in a higher temperature range, and no decrease in it at higher 
temperatures It is not clear if those effects are the results of nonlinear 
phenomena or of harmonic vibrations [0, A calculation of transport in 
a vibration model is a good way of testing whether the main mechanism of 
heat conductance in glasses is obtained from the linear harmonic vibration 
modes, or if nonlinear scattering processes have to be invoked. 

Since vibration Hamiltonians are elastically stable there are special con- 
straints on the dynamics. Therefore, though the dynamical equations are 
similar to the Anderson problem of a quantum particle in a random poten- 
tial there are some major diffrences due to these constraints. This model is 
a random Laplacian model, with positive off diagonal coefficients. It leads to 
the existence of a global translation mode at zero frequency. This features 
has a decisive influence on the entire spectrum of vibrations at low frequen- 
cies. Other distinguishing features are the vector character of the vibrations, 
and in the case of a two (or more) component systems, random masses. We 
will comment in more detail on these differences in section 2. 

The focus of this study is a scalar lattice vibration model with random 
elastic coefficients, and equal masses. This model disregards some of the 
important features of glass vibrations. Nevertheless, in a previous work which 
focused on the eigenstates , it was shown that this model contains many of 
the interesting features of glass vibrations, and is qualitatively different from 
the usual electronic models. The present work reveals the phase diagram and 
the conductance properties of the disordered lattice model. We note that the 
model under discussion also describes a free quantum particle with disorder 
in its transport behavior. 

In our study we find several new qualitative features that are absent 
in the classical Anderson problem. Some theoretical studies were vibration 
problems with disorder, such as P, |^, have concluded that they belong to 
the same universality class as the Anderson model. This is not the case for 
in the model studied here and in 

We outline the main differences: In the three dimensional electronic An- 
derson model there are two phases: diffusive with extended states and lo- 
calized. In the vibrational case there are four phases. In addition to the 
two phases familiar from the Anderson problem there is an insulating sub- 
exponentially localized phase with multifractal normal modes, and a strongly 
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conducting phase, with states characterized by weak scattering, and a well 
defined wavelength. The phase diagram is displayed in fig. 0, and described 
in detail in section ^. 

In two dimensions the Anderson model is insulating with exponential lo- 
calization and all the states are multifractal, while in the vibrational case 
exponential localization is also observed, but a transition occurs between 
multifractal and extended modes, accompanied by a singularity in the local- 
ization length. This is discussed in detail in our paper on the two dimensional 
case 0. 

The former study |^ concentrated on a numerical study of the properties 
of the normal modes. In the present paper we apply the transfer matrix 
method ||^ to the vibration model. This method has the advantage of point- 
ing out precisely the transition between conducting and insulating behavior 
through localization. The finite size scaling analysis which has been used 
earlier to calculate the physical quantities now reveals the presence of addi- 
tional length scales. The method also provides a direct tool for calculating 
the conductance of finite systems. The various numerical methods which 
are used to analyze the Anderson model, and the results obtained have been 
recently reviewed by . 

The structure of the paper is as follows: in section |^ we define the model, 
the numerical method, and explain how the numerical data are used to study 
the system. Section ^ describes in detail the results in the three dimensional 
case, emphasizing the different qualitative behaviors. In section ^ we apply 
the results of section |^ to calculate the thermal conductivity as a function of 
temperature. The paper is concluded with a discussion of the relevance of 
our results to glasses. 



2 The disordered harmonic model and elastic 
energy transport 

The vibration problem in glasses is a vector problem which is defined for a 
disordered lattice where both structural and elastic disorder play a major 
part. In this paper we treat a simpler variant of this problem: the scalar vi- 
bration model. The model under consideration is defined on a (i-dimensional 
simple cubic lattice with periodic boundary conditions. It consists of a scalar 
field which is subject to elastic equations of motion. For a normal mode of 
frequency the field equations read 

-u^{x)= k.^Am-m]. (1) 
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gions of different quahtative transport behavior, as a function of the fre- 
quency u, for a fixed disorder distribution. The small u region below the 
dashed line is present only in finite systems. The descriptions given in italics 
refer to normal modes, analyzed in [^. The symbols appearing in the dia- 
gram are the finite size localization length C,w for a system of width W, the 
conductance Gw of a finite cubic system of side W, and the conductivity g. 
They are explained in detail in section ^. 
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where x and y denote lattice sites. The fc's are symmetric (A;^,^' = 
positive, independent and identically distributed random numbers. A similar 
model was introduced recently by [TT| to describe boson peaks in glasses. 

This model is in fact a realization of the Anderson tight binding model 
with special constraints. The diagonal terms are the sum of the off diagonal 
terms in the same row, and the off diagonal terms are real and positive 
number. The same model describes free particle. 

The focus of the present study is the transport properties of this model. 
For this purpose it is convenient to single out a preferred direction, and label 
the coordinate in this direction by t. We define a column vector Vt whose 
components are the values of (p at the sites on the hyperplane labeled by t. 
It is then possible to use eq. (|I]) to express Vt+i in terms of Vt and Vt-i. This 
linear relation defines a random transfer matrix T by 

Vt+l \ _rj, f Vt 



Tti ' . (2) 

Vt J \Vt-lJ ^ ' 

The transfer matrices T have the same non-zero elements as those of the 
Anderson tight-binding model. 

Following the standard methods developed in the context of disordered 
electronic systems, we study the conductance properties using the Lyapunov 
spectrum of the product 

••■TiTi_i---Ti . (3) 

Due to the left-right symmetry of the problem, the Lyapunov exponents come 
in pairs of opposite signs. For any positive value of a;^, and for no n- vanishing 
disorder all the exponents are non-zero, signaling that in a quasi-lD geometry 
the system is always insulating. The smallest of the non-zero exponents may 
be identified with the best conducting channel, so that it gives the inverse of 
the finite size localization length where W is the transverse size of the 
system. 

We study as a function of W , with the aim of finding the asymptotic 
behavior for large W . The following types of asymptotics are known to occur 
in models of electronic localization 

1. The localization length saturates to a finite value ^ ^oo = const . 
This is the localized phase which is expected to characterize systems 
with strong disorder, and low dimensionality. 

2. The critical state, E,w ~ W . This behavior is commonly associated 
with the delocalization transition, using the scaling hypothesis. 
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3. The diffusive phase, ~ DW^ ^ This phase is expected to have 
Ohmic behavior, 

G ~ gW^-^ , (4) 

and the conductivity g is proportional to D, 

Since the transverse size of the system which can be studied numerically is 
limited, the extrapolation to infinite size can be improved significantly using 
finite-size scaling (FSS) 0, which assumes a scaling form 

= UfiW/U , (5) 

with different scaling functions above and below the transition. The local- 
ization length is found numerically to diverge as 

^oo ~ |0- - (Tel'" , (6) 

in three dimensions and as 

~ exp (a--^') , (7) 



in two dimensions where a measures the disorder strength, and dc is the 
critical value of disorder in three dimensions. 

The average conductance G is given, using the the entire set of Lyapunov 
exponents •jj, by 

^ = ^cosh^(7,iy) ' 

for a hypercubic sample [|1^. We will show below that in the diffusive state, 
since the spacing between the Lyapunov exponents is slowly varying, the 
conductivity g is indeed proportional to D. However, this is not always 
valid, and we show below that the relation between g and D breaks down 
when scattering is too weak. 

Here, we study the transport properties for different values of u, rather 
than disorder. As expected, we find that small u corresponds in many ways 
to weak disorder; the basic reason for this is that the uncorrelated disorder 
is averaged out in its effect on large wavelength modes, which correspond to 
small u. However, the mapping between small uj, and weak noise is exact 



only in one dimension We find that the behavior as a function of u is 
considerably more involved than that which is known for electronic models, 
and which has been reviewed above. We find several new phases, and con- 
clude that in three dimensions the simple scaling picture of the localization 
transition does hold for our model. As a result, FSS does not hold in the 
simple form (|^), and has to be modified. The present study confirms and 
adds further information to the new features which were discovered in the 
geometry of the eigenmodes of the same model |^ . 



6 



3 Transport properties in three dimensions 



As described is section ^, the various phases of the system can be character- 
ized by the asymptotic behavior for large transverse width W of the the finite 
size locahzation length which is the inverse of the smallest non-negative 
Lyapunov exponent and of the conductance G which is defined in terms of 
the Lyapunov exponents in The results of this section were obtained 
for the model (|l|) in a bar shaped geometry, with a square cross section of 
width W, and a fixed distribution of the values of the elastic constants kg^^j: 
A binary distribution with Prob(/c = 0.1) = 0.8 and Prob(/c = 1) = 0.2. We 
computed the Lyapunov exponents of the product (0) for different values of 
W, 5 < W < 16. The number of realizations which were multiplied was 
between 4 x 10^ for the smallest widths, and 2.5 x 10^ for the largest. 

We studied the transport properties for varying u. These results are 
summarized in the phase diagram displayed in fig. |l|. We sketch the major 
features of the phase diagram. For frequencies which are higher then ui 
we observe a normal localized phase with exponential localization. However 
in the transition range near Ui, the usual finite size scaling ansatz fails, 
and a modified finite size scaling with two parameters is needed to rescale 
the data. For smaller frequencies in the range uj2 < uj < uJi we observe 
a divergence of the finite localization length as a non-trivial power of the 
system size in the entire range available to us. This phase is interpreted as 
being sub-exponentially localized. In the range < u < 002 we observe 
the normal diffusive behavior familiar from electronic models, but finite size 
scaling has again to be modified near the transition. These effects are related 
to a continuous change in the multifractal exponents of the normal modes in 
this range, instead of a sharp transition. 

At frequencies which are lower than ^3, the normal diffusive behavior is 
replaced by an irregular dependence of the transport properties on the system 
parameters. The diffusion coefficient becomes hard to define numerically, and 
a variety of types of size dependence of the conductance is observed. Since 
uj < UJ3 is the range corresponding to weak scattering as found in [^, these 
effects can be interpreted as residuals of the ballistic propagation between 
scattering events. We expect a long crossover region before an ultimately 
diffusive behavior can be observed. 

The rest of this section is devoted to a detailed description of the phases 
in the order they appear with decreasing u. 
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3.1 The normal localized phase 

The system is localized for values of u;^ > 0;^ = 1.1. This phase is character- 
ized by the existence of a finite limiting localization length 

^00 = „lim <oo , (9) 

which we conventionally identify with the range of the Green's function of 
eq. (0) 

eoo = ,Jim --^{\og\g{x,y;iu)\) . (10) 

\x-y\^oo \X — y\ 

We found that in the localized phase the dependence of on 1/W is 
well approximated as linear. This function can then be linearly extrapolated 
to 1/W = 0, giving an estimate for l/.^oo- The values of ^00 thus obtained 
are listed in table for various values of tu^. 

As can be seen from table |1|, the value ,^00 increases as u"^ decreases as 
it starts to diverge near the transition. However, because of the limited 
size of the systems which were measured, near the transition the precision 
of the measured ^00 deteriorates, until all precision is lost. As discussed in 
section |^ above, this difficulty is commonly dealt with by invoking the FSS 
hypothesis (|^), which can be written equivalently as \og{^w/W) = g(\ogW — 
log ^00); Coo(i^) is then determined so that the various curves of measured 
log{C,w /W) versus logVT collapse. 

This procedure works well for u"^ > 1.14. However, for uf = 1.1 < u"^ < 
1.14 the numerical curves do not collapse implying that FSS is violated. It is 
nevertheless possible to describe the measured values of for all u"^ using 
a modified form of FSS 



K^')/ ( 77^ h (11) 



the values of log a and log b thus obtained are also reported in table |1|, and 
the resulting data collapse is displayed in fig. ^. Unfortunately the fitting 
parameters a and b are not directly related to the localization length, and 
therefore the modified FSS is not useful for a measurement of ^00 near the 
transition. The physical interpretation of a and b remains unclear. For 
00"^ > 1.14, where b = 1, a may be again identified with ^oo, which gives 
values in accord with those obtained from direct extrapolation. 

As discussed above, ^00 diverges at the transition, but it is difficult to 
measure this divergence accurately. Nevertheless, the values of ,^00 obtained 
by direct extrapolation are well in accord with a power like divergence 

^r^iu-U,)--^ (12) 



8 





c 

Soo 


log 


10 ^ 


logio ^ 


2.0 


0.377 








1.6 


0.131 








1.4 


0.226 


-0 


16 




1.3 


0.259 


-0 


21 




1.25 


0.160 


-0 


02 




1.2 


0.094 





22 




1.17 


0.073 





33 




1.15 


0.047 





56 




1.14 


0.033 





73 




1.135 


0.026 





85 




1.13 


0.019 





97 


-0.003 


1.125 


0.013 


1 


11 


-0.007 
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Table 1: The values of the localization length ,^oo as determined by extrap- 
olation to ly = oo, and the parameters a and h of the modified FSS (H), as 
a function of cj^ in the localized phase in three dimensions. The values miss- 
ing from the table could not be determined significantly numerically, except 
when a value is given for a and not for 6, where h was fit to 1. The different 
symbols were obtained for different values of a;, all of which are given in 
table |1|. 
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Figure 2: Modified finite size scaling in the localized phase in three dimen- 
sions. W is the width, and C,w is the finite size localization length; the 
logarithms are in base 10. The data displayed are for 1.105 < < 1.4. The 
FSS parameters a and b as functions of uj are reported in table |1] 
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with ui ~ 1.7. This value is within the numerical errors in agreement with the 
value of z/ = 1.5 reported for the three-dimensional Anderson transition ||13 |. 
On the other hand, the fitting parameters a and b increase very slowly near 
Ui, and their rate of divergence was impossible to determine to a reasonable 
accuracy, from the numerical measurement. The breakdown of FSS implies 
there is at least one more length scale, in the system in addition to ^oo(and 
W); unlike ^oo, ^* does not diverge at the transition. 



3.2 The anomalous localized phase 

This subsection describes the properties of the system for = 0.96 < tu^ < 
ujf = 1.1. In electronic models the Anderson transition in three dimensions 
is the boundary between localized and conducting diffusive phases. In the 
present model the qualitative behavior is markedly different. One surprising 
feature is that the finite size localization length for a given W viewed as 
a function of u"^ has a maximum at the the transition point u"^ = 1.1. That 
is, the conductivity of a finite sample becomes worse when u is decreased 
from the transition, the opposite of what would occur if the transition were 
into a conducting phase. This trend continues until a minimum is reached 
at tu^ = 1.0, and the usual situation whereby increases with decreasing 
is restored (see fig. 

The behavior for a fixed as a function of W is also unusual. C,w does not 
saturate to a finite value, and continues to grow with W as some non-trivial 
power 

^w^-W^ , 0</3<l . (13) 

Some examples are given in fig. ^, where a clear power law behavior is ob- 
served, in contrast with the other phases. 

The measured values of (3 are reported in table |^, and the minimum occurs 
at cij^ = 1. An interesting feature is that the at the lower bound of this phase 
(3 is still significantly smaller than 1; we discuss below the implication of this 
fact on the nature of the second phase transition into the diffusive phase. 
The asymptotics of C,w given in eq. (|1^) shows that the Green's function in 
this phase is certainly not localized with a finite localization length ^oo in 
the sense defined in (pUj). However, it is compatible with a weaker 'stretched 
exponential' localization 



{\Q{x,y)\) ~ exp 



\x-y\ 



(14) 



with < a < 1. If the three dimensional system indeed behaves as in (li) 
then we can expect the same behavior in the bar geometry which we stud- 
ied for distances \x — y\ which are small compared to ^w, which marks the 
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Figure 3: Finite size localization length divided by the width W ^ as a 
function of the frequency squared tu^, for widths = 5, 8, 12, from top to 
bottom. The vertical lines show the limits of the anomalous localized phase. 
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Figure 4: The base 10 logarithm of the Finite size locahzation length 
divided by the width W, versus log^Q W, shown for three frequencies in the 
anomalous localized range. The full, dashed, and dot-dashed lines refer to 
uj^ = 1.08, 1.00, 0.97 respectively. 
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Table 2: The localization exponent (3 and the anomalous localization length 
as a function of cj^, in the anomalous localized phase. 



crossover to true exponential decay. liW^ we can estimate by de- 
manding that the two types of decay are of the same order of magnitude for 
distances of order W , that is, 



exp 




exp 



W_ 



(15) 



This yields the estimate 

which agrees with ( p!4|) identifying P = 1 — a. 

The delocalization at uj ^ uJi from below, as described in this section, 
occurs because /3 — > 1, and the length remains finite. It seems natural 
therefore to identify C,* with the additional length scale present in the normal 
localized phase. 

The transition in U2 is unique since a{uj2) > 0. One has a stretched 
exponential correlation with an internal length C,* which should diverge to 
infinity at the transition. This is analog to the usual process of delocalization. 
The value of ^* can be calculated from the numerical results using eq. (p!5D, 
and the values obtained are reported in table ||. However, the measured 
values of are inaccurate, and should be considered only as indicative. If 
these values are trusted, the conclusion is that the divergence of ^* near 002 
is very slow, but a quantitative statement is impossible to make. 

In the analysis of the normal modes in parameter range described in 
this section could not be observed with a good enough resolution. A more 
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detailed study enabled us to observe this range with better statistics, and the 
conclusion is that the states are multifractal in this range, with exponents 
changing continuously, in a way similar to the change of the exponent (3 of 
the divergence of the correlation length. These results and the relationship 
between them and the present study will be given in a separate paper . A 
possible way to describe the phase described in this section, is that it is an 
intermediate regime, where the change in the fractal exponents smears the 
Anderson transition. 



3.3 The diffusive phase 

This is a conducting phase with Ohmic behavior, observed in the parameter 
range a;| = 0.45 < u;^ < = 0.96. The characterizing feature of this phase 
is the existence of a well-defined nonzero limit 

lim ^ = D . (17) 

D is approximately proportional to the conductivity of the system, as is 
discussed in detail below. The lower limit of this phase tus is not very well 
defined, because finite size effects are stronger for small u. 

As in the normal localized phase (see subsection |3.1| ) there are two meth- 
ods in principle to calculate D from the raw numerical data. One may 
regard D^r = ^w/W^ as a function of 1/W and extrapolate to 1/W = 0, 
which should work well far from the transition point u)2, where D is not too 
small, or, near the transition D may be inferred from FSS, eq. (|^), where in 



this case the scaling function / grows quadratically for large arguments |jT3 
However, as in the normal localized phase, the standard form of FSS breaks 
down near the transition and has to be replaced by the modified FSS (pJ^, 
which is not useful for the determination of D. Again, this breakdown implies 
the presence of additional length scales. 

Since FSS is of limited use, the measured values of D reported in table | 
were obtained by straightforward extrapolation. The values approach zero 
when u; ^ co'2 as in the usual Anderson transition. The decay rate can be fit 
quite convincingly with power law 

D{lj) ~ (CJ2 - ^^T^ , (18) 

with z/3 = 1.5 ± 0.1, again in agreement with the Anderson transition, 
within numerical error with Ui which measures the divergence of the local- 
ization length near the boundary between the normal localized phase, and 
the anomalous localized phase [see eq. (P!^)]. 



15 





D 


logio a 


logio b 


9 


0.50 


0.891 


0.786 






0.60 


0.657 


0.00 




0.554 


0.70 


0.423 


0.20 




0.354 


0.75 


0.283 


0.35 




0.245 


0.80 


0.173 


0.57 




0.148 


0.82 


0.127 


0.68 




0.115 


0.85 


0.076 


0.90 




0.070 


0.87 


0.050 


1.08 




0.047 


0.89 


0.032 


1.29 




0.030 


0.90 


0.025 


1.40 




0.024 


0.91 


0.019 


1.55 




0.018 


0.92 


0.014 


1.74 




0.012 


0.93 


0.008 


2.03 




0.005 


0.933 


0.006 


2.10 


0.005 


0.004 


0.935 




2.18 


0.008 




0.938 




2.28 


0.014 




0.94 




2.36 


0.019 




0.943 




2.49 


0.028 




0.945 




2.60 


0.036 




0.946 




2.68 


0.041 




0.947 




2.71 


0.045 




0.95 




2.90 


0.063 




0.953 




3.14 


0.088 
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Table 3: Same as table |l| for the three-dimensional diffusive phase, except 
that D is given instead of S,ooi and the conductivity g is also given. It is 
evident that D and g are close throughout this phase. 
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Figure 5: Same as fig. ^ for tlie diffusive phase in three dimensions. The 
frequency range is 0.5 < a;^ < 0.955. a and b as functions of uj are reported 
in table ^ 

The values of log a and log 6 obtained from the modified FSS ( pA] ) are 
also reported in table ^, and the data collapse obtained after scaling the raw 
data is shown in fig. ^ For values of uj where b = 1, the parameter a may 
be identified with ^^o = with good agreement with the values of D 

measured by direct extrapolation. 

As can be seen from table ^ the parameters a and b increase in absolute 
value as u is increased toward uj2, and diverge slowly. This implies an in- 
teresting behavior of the localization properties near uj2'- Suppose that the 
asymptotic behavior of the scaling function / in ([llD for small argument is 
power hke, 

f{x) ~ , for a; ^ . (19) 
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Taking a value of a; close to a;2, so that a is very large, one finds from eqs. (|Tl|) 
and (|19D that 

~ (20) 

for W a. This behavior has to match the one observed at the lower edge 
of the anomalous localized phase, which implies the identity 

p = p(uj2) ~ 0.88 . (21) 



The dependence on W shown in (|20|) means that very close to the transition 
the system looks localized for small W because ^w/W is decreasing, before 
the crossover to the asymptotic linear growth of C,w/W. This scenario is 
corroborated by the minimum in the curve shown in fig. ^ 

Since this is a conducting phase, it is instructive to compare the pa- 
rameter D with the conductivity g. The calculation of the latter from the 
Lyapunov spectrum may be achieved using eqs. (H) and (^). We remark that 
the dimensionless conductance G calculated from (H) is sometimes called the 
transmission coefficient, since it is obtained by regarding the system as placed 
between ideal leads and providing W^'^ channels which conduct energy be- 
tween the leads. The precise dependence of G on the channel transmissions 



is model dependent, and (||) refiects one possible choice |jT4 . 

The reason that D gives a good estimate for g in the diffusive phase is 
that the Lyapunov spectrum in this phase vanishes linearly at zero [0 , and 



the spacing between the Lyapunov exponents changes slowly, as exemplified 
in fig. 1^. This is in contrast with the Lyapunov spectrum in the low frequency 
range discussed below. 

To see how g and D are related let us assume that the W'^ Lyapunov 
exponents are evenly spaced, with spacing l/C,w- The formula for the con- 
ductivity (H) involves a sharp cutoff at 7„ ~ W, so it basically counts the 
number of exponents smaller than 1/W. Therefore we get G ~ ^wf^- 
Since ~ DW^ and G = gW, g and D are proportional. In table |^ we give 
the measured values of g. The agreement between D and g remains good 
throughout the diffusive phase. 



3.4 The strongly conducting phase 

In this subsection we describe the properties of the system for uj'^ < = 0.45. 

We find that for very small u short range disorder becomes irrelevant, and 
the vibration modes behave very similarly to modes of a pure system with a 
well-defined wave vector k and an acoustic dispersion relation u = c\k\. Such 
modes exist even in one-dimensional systems, but are restricted to a band of 
width W~^^'^ around a; = 0. Similarly in three dimensions, although for any 
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Figure 6: A typical Lyapunov spectrum in the diffusive phase, with the 
exponents decreasing hnearly to zero. 7^, i — 1, • • • , 144 are the positive 
Lyapunov exponents for a bar of width 12, at cu^ — 0.6. 
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positive uj the transport is ultimately diffusive, there is a band of ballistic 
modes with a well defined wave vector near = 0, which disappears in the 
thermodynamic limit. The existence of this band may be simply understood 
as a consequence of the divergence of the mean free path £ as a; — >■ 0. For any 
finite system there is a positive frequency where I becomes larger than the 
system size, and for smaller frequency a wave propagates through the system 
with essentially no scattering. 

Outside the ballistic band, and for uj < lo^-, the normal mode analysis 
of shows that the modes experiences scattering which destroys the unique 
wave vector dependence. However, the absolute scale of the wave vector, the 
wavelength^ is still well defined, and it is significantly larger than the lattice 
spacing. The wavelength is defined up to some width, which is commonly 
associated with the inverse of the mean free path. 

The transport analysis of the present study reveals that a qualitative 
change in the structure of the Lyapunov spectrum occurs at ujj,. In fig. ^ 
a typical Lyapunov spectrum of the low frequency phase is displayed, to be 
compared with fig. ^. The Lyapunov spectrum at small lo is characterized 
by the presence of large gaps; these obviously rule out the possibility of a 
description of the conductivity using a single parameter. The conductivity 
of the system depends sensitively on the system size, as well as the other 
parameters of the problem. 

The origin of the gaps in the Lyapunov spectrum can be traced to the 
behavior of a pure system which is extended infinitely in the longitudinal 
direction. In such a system, for every longitudinal wave number fcy, there are 
W"^ modes with 

^l = cl{k\^kl) (22) 

where /cn, 1 < < are the possible values of the transverse wave vector. 
These represent W'^ channels, which are perfectly conducting in the pure 
system. However, for a given lo one will have less then W'^ solutions since 
some of channels may be blocked, when the transverse vector is too large 
(transverse confinement). Therefore the conductivity of the pure system 
depends discontinuously on uj: When uj is increased beyond the threshold 
values the conductivity jumps by an integer amount. 

This picture describes quite accurately the conductance G for small fre- 
quencies ijj <^ cufe. As ijj is increased the Lyapunov spectrum is less well 
described by the one of an effective pure system, but gaps and irregulari- 
ties in the spectrum persist up to ws, reflecting a no n- vanishing transverse 
confinement. 

Consider the conductance G as a function of lo for a fixed system size W . 
As just explained, for small lo system size the conductance is ballistic, and the 
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Figure 7: A typical Lyapunov spectrum in the strongly conducting, for a 
bar of width 12, at iJ^ = 0.2 (compare with fig. p.) 
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value is proportional to the number of open channels, which is proportional 
to {ujWy. This behavior persists as long as W is smaller than the mean free 
path i. Since the mean free path diverges as £ ~ uj~^ in three dimensions 
P, |ll|, the crossover occurs at cui, ~ therefore G{ujb) ~ W^/\ 

Since the system is expected to be diffusive for any positive u, the de- 
pendence of G on should change from the quadratic dependence of the 
ballistic range to an ultimately linear one, G ~ gW. If we assume a simple 
crossover from ballistic to diffusive behavior of G when W ^ i, continuity of 
G inW implies that {uJbWy ~ W^^^ ~ g{uJb)W, and the conductivity would 
diverge as g^uj) ~ uj~^. However, the conductivity should diverge at the same 
rate as the mean free path, i.e., as It follows then that there is a range 
in W where G grows stronger than linearly but weaker than quadratically, 
before the final diffusive behavior sets in. 

These different behaviors and crossovers are represented in the conduc- 
tance as measured numerically, and presented in fig. ^. The conductance at 
small u! is very well described by that of an effective pure system with sound 
speed Cs = 0.35, which is between the two possible values for the elastic 
constant, 0.1 and 1 (recall that the masses were set to 1). The agreement is 
quantitative as is demonstrated in fig. ^, where the measured conductivity is 
shown against u'^W'^, and compared with that of the effective pure system. 
At the frequency range ujf, < uj < the functional form of G is harder to 
describe. Gw indeed grows faster than linearly and slower than quadrati- 
cally, but as a result of the transverse effects discussed above, the growth 
rate is not monotonic and subject to sharp changes. Thus, we were unable 
to measure precisely the growth rate of G as a function of W for a fixed u. 
On the other hand it has been possible to measure the growth rate of the 
maximal value of G for a fixed W 

maxGv^M ~ Vr^'S (23) 

which is the smallest possible growth rate, since ma.x^ Gw{^) > G{ujb) ~ 
W"^!"^ . The ultimately linear growth of G as a function of W could not be 
observed in our numerics in this range for our system sizes. 

4 Thermal conductivity in three dimensions 

The results of the previous section can be used to calculate the overall heat 
conductivity of a three-dimensional disordered harmonic lattice, in the pres- 
ence of a temperature gradient. The heat flow through a sample connected 
to two infinite ideal leads at different temperatures has been recently de- 
rived [0. From this, the thermal conductivity k may be calculated, taking 
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Figure 8: The conductance G as function of the frequency squared 
uj"^ for different bar widths. From bottom to top the widths are W = 
6, 8, 10, 12, 14, 16. The dependence on W for uj"^ > 0.45 is hnear. 
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Figure 9: The same data as in fig. |^ but shown as a function as u'^W'^. The 
top curve is the result that would be obtained for a pure system with sound 
speed Cs = 0.35. 
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the temperature difference to be small, giving 



fl 1 /-oo 

= I, duJudTm{^)G{uj) , (24) 

where G{uo) is the frequency dependent transmission rate, discussed in the 
previous section, and riT{ijj) = l/(e^'^/"^ — 1). Below we use units such that 
Planck's constant and Boltzmann's constant are both equal to one. 



The distribution function ?7t(i^) which appears in (pif) diverges like T/uj 
for small u, and has an exponential cutoff at a; ~ T. Hence for the purpose of 
a qualitative estimate we can replace rj^iuj) by T/u, cutting off the frequency 
integration at = CT, with C an order 1 constant. This approximation gives 

K ~ — / dujG(Lj) . (25) 
W Jo 

Thus, as the temperature increases, modes with higher uj contribute to k; if 
CT is in the frequency range of one of the phases, then some of the modes of 
this phase contribute to k, the phases of lower uj contribute fully, and those 
of higher uj do not contribute at all. 

It follows that when T <ti uj^ only the strongly conducting phase should 
be taken into account. If the system is small enough so that T ^ uj^, the 
conductance is carried by ballistic modes, and eq. (pSf) gives 



K ~ WT^ , (26) 

as in a perfect lattice. When T > this form is no longer valid, and 
scattering should be taken into account. Unfortunately, the numerical results 
provide limited information about the asymptotic behavior of G in this range, 
and we are forced to make an ansatz about the form of G which takes into 
account the numerical results and theoretical predictions, such as 

uj^W^ < < cjfo ~ W-^/^ 
G ~ { iy3/2 < ^ < ^ _ (27) 

UJ~^W UJc < UJ < UJs 

This ansatz gives a linear dependence of k on T in the intermediate frequency 
range, with a coefficient proportional to W^^^ and a fast saturation at higher 
temperatures. 

At temperatures greater or of the order of uj^ the contribution to k, of the 
strongly conducting phase is therefore anomalous growing with the system 
size as W^^^. The exponent of the anomaly should be considered as only in- 
dicative, but the contribution of the ballistic band itself already gives a W^^^ 
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anomaly. The divergence of the heat conductivity of the harmonic modes 
at low frequency is quite robust, and has been observed also in molecular 
simulations of glasses [Q. Nonlinear interaction are invoked to obtain heat 
conductivity which remains finite in the thermodynamic limit. This can be 
done phenomenologically by replacing the system size dependence of k by 
the inelastic mean free path. 

In contrast, the contribution of the diffusive phase is quite simple: 

^diffusive ~ / D{uj) , (28) 

which gives Ohmic, system size independent, conductivity, increasing linearly 
at first, and saturating as T ~ uj2- This contribution, however, becomes 
negligible when the system size is taken to infinity. Finally the contribution 
of the localized phases to k decays exponentially in the system size, and these 
phases are insulating. 

We also used eq. to calculate k(T, W) directly from the numerical 
values of G{ijj)] the results are displayed in fig. |T0|, scaled with W^^"^. The 
integration over uj suppresses the finite size effects, and it can be seen that 
the W^^"^ scaling slightly overestimates the measured anomaly. As expected, 
one observes in fig. llO] a increase for small T, followed by a linear increase 



for intermediate T and a saturation. A power law fit for the saturation yields 
an W^^^ anomaly, but we regard this value as indicative only. 



5 Conclusions 

We wish to comment on two central issues related to our model. First, in 
the scalar vibration model the delocalization transition is smoothed in the 
sense that there are two critical points, and there is a continuous change 
of the multifractal exponents near the transition. Moreover, far from the 
transitions one can observe the classical finite size scaling behavior similar 
to the Anderson scaling. 

Our second principal conclusion is that the heat transport of modes of 
low frequency has some peculiarities probably related to the existence of a 
well defined wavelength. The low-temperature conductivity due to this band 
reproduces some of the experimental phenomena in glasses, such as a linear 
dependence on the temperature, and an eventual saturation (plateau). In 
recent papers [^, it has been argued that the heat transport in silicon glass 
is governed by harmonic modes; our model reproduces quite convincingly 
several of their results. In particular, there is no need to invoke non-linear 
processes such as two level systems to explain the temperature dependence 
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Figure 10: The thermal conductivity k of a cube of side L divided by the 
square root of L, as a function of the temperature T, shown for different Ls. 
The values of L shown are (from to to bottom) 6, 8, 10, and 12. 
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of the heat conductivity. In all cases, however, the harmonic approximation 
yields a conductivity which diverges as some small power of the system size, 
and this feature must be corrected by the nonlinearity of the interaction. 

The disordered lattice model is thus found to be a useful tool to study 
heat transport of glasses, and it will be very interesting to find out what will 
happen with a vector vibration field and with random masses. 
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